Nonlinear magnetoplasmons in strongly coupled Yukawa plasmas 
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The existence of plasma oscillations at multiples of the magnetoplasmon frequency in a strongly 
coupled two-dimensional magnetized Yukawa plasma is reported, based on extensive molecular dy- 
namics simulations. These modes are the analogues of Bernstein modes which are renormalized by 
strong interparticle correlations. Their properties are theoretically explained by a dielectric function 
incorporating the combined effect of a magnetic field, strong correlations and finite temperature. 
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Two-dimensional (2D) systems with strong correla- 
tions, i.e. with the interaction energy exceeding the ther- 
mal energy, show a number of unusual properties ranging 
from anomalous phase transitions (Kosterlitz-Thouless 
scenario), to anomalous transport, e.g. [TJ [2]. Exam- 
ples include electrons on helium droplets and in quantum 
dots [3], ions in traps [4] and dusty plasmas, e.g. [5[|6]. 
A particular correlation effect is observed in the collec- 
tive oscillation spectrum of 2D Coulomb and Yukawa liq- 
uids [5j [7J |8] where a transverse shear mode has been pre- 
dicted by Kalman and Golden [9l [10] . This mode does 
not exist in an ideal system but has experimentally been 
observed in strongly coupled systems [TTJ [12] . Similarly, 
in the presence of a strong magnetic field a 2D plasma 
shows two coupled modes - a magnetoplasmon and a ma- 
gentoshear mode (upper and lower hybrid modes) which 
have recently been studied in detail in the liquid and 
crystal phases [BHH]. 

Here we show that a strongly correlated 2D Yukawa 
plasma (2DYP) possesses additional collective modes. 
We demonstrate that these are related to Bernstein 
modes (BM) [13] which are well known in high- 
temperature classical, e.g. p~4j[15] and quantum plasmas, 
e.g. [l6j[T7]- However, these are all nearly ideal plasmas. 
Here we demonstrate that such modes can also exist in 
strongly correlated systems but they have a number of 
peculiar properties arising from the combined effect of 
correlations and the magnetic field. 

Model and simulation method. Our 2DYP con- 
sists of N particles in a quadratic monolayer (with pe- 
riodic boundary conditions), subject to a magnetic field 
perpendicular to the plane and interacting by Yukawa 
forces. The coupled Newton's equations are (i = 1 . . . N) 

m r i = |r iX B-Q 2 ^(v^^ , (1) 

r=\fi-fj\ 

with the Debye screening length A^ and charge Q. 
Here we define the Wigner-Seitz radius a = (no7r) -1 / 2 
(where no is the areal density) and the frequency 



ujo = (27rQ 2 no/ma) 1 / 2 , to be distinguished from the k- 
dependent 2D analogon of the plasma frequency, 0J p (k) = 
[27rQ 2 nok 2 / (m\fk 2 + ft 2 )] 1 / 2 . In thermodynamic equilib- 
rium, the system is fully characterized by three param- 
eters - the dimensionless inverse screening length ft = 
a I Ad, the Coulomb coupling parameter T = Q 2 / \aksT) 
and the dimensionless magnetic field strength, (3 = 
uj c /ojq, with the cyclotron frequency uj c = QB/mc. Eq. 
([!]) is solved by a standard microcanonical MD technique, 
with data acquisition starting after the system has equi- 
libriated to the desired temperature (T). The longitudi- 
nal and transverse collective excitation spectra, L(k,u) 
and T(k,(j), respectively, are computed from the micro- 
scopic longitudinal and transverse current fluctuations, 
see e.g. [18]. Here we use N = 4080 particles. 

Numerical results. A first overview of the simula- 
tion results is given in Fig. [I] showing the evolution of 
the collective oscillation spectrum with the coupling pa- 
rameter T at a fixed field strength, /3 = 0.5. At low wave 
numbers k one clearly observes two peaks, one of which 
starts at uj = uj c = 0.5 and a second at very low fre- 
quency. When T is increased, cf. Fig. [T^l-c, the peaks 
become more pronounced and extend to higher k. These 
are the well-known magnetoplasmon (MP) and magne- 
toshear modes - the latter being a pure correlation effect 
as noted above [9]. When the coupling is increased be- 
yond the crystallization point (when the system trans- 
forms into a microcrystalline structure in our simula- 
tion), cf. Fig. [l]d, the modes transform into the mag- 
netophonon spectrum, e.g. [6j[8]. The properties of these 
two modes are well described within the quasilocalized 
charge approximation (QLCA) and harmonic lattice the- 
ory and very well agree with MD simulations; thus it 
appears that the collective excitations of strongly corre- 
lated Coulomb and Yukawa systems are fully understood. 

However, this is not the case. Let us look at frequencies 
above the MP. For T = 80, Fig. [T]b, a third maximum 
appears in the spectrum around 1.3cjo, starting at a wave 
number of about ka = 1.5. This peak is very broad and 
strongly overlaps with the MP. An analysis of the spec- 
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FIG. 1: (Color) Collective excitation spectrum, L(k,cj) + 
T(k,cj), of a Yukawa plasma with k — 2 and /3 = 0.5, for 
four values of the coupling parameter given in the figure. The 
amplitudes of the spectra are normalized by their peak values. 
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FIG. 2: (Color) a, c, d: Longitudinal magneto-oscillation 
spectrum L(k,uj) of a Yukawa plasma for various field 
strengths and coupling parameters, b: Transverse spectrum 
T(k,uj) for the parameters of part a. The lines are vertically 
offset for better visibility. 



trum for T = 200 and 1000 reveals that this peak exists 
there as well - it is now more clearly separated from the 
MP, cf. Fig. [l]3,d. To explore the origin of this peak we 
have performed a series of simulations for a broad range 
of field strengths, coupling parameters and screening pa- 
rameters. All the results confirm the appearance of not 
just one peak but of a whole series of additional collec- 
tive modes. Some representative results are collected in 
Figs. [2] and [3) Both fi gures clearly show up to four addi- 
tional plasmon modes, which are equally spaced, where 
the spacing increases with B. 

These additional peaks and their monotonic depen- 
dence on the cyclotron frequency strongly resemble Bern- 
stein modes [T3] . Furthermore, the conincidence of lon- 
gitudinal and transverse current fluctuation spectra, cf. 
Fig.[2j top row, clearly confirms the circular polarization. 
In fact, electrostatic modes with frequencies around mul- 
tiples of uj c are well known in many ideal classical and 
quantum plasmas, see above. But here they are observed, 
for the first time, in a strongly coupled system. Note that 
neither cold fluid theory nor the QLCA have revealed 
such modes because they neglect finite temperature ef- 
fects, which are crucial for the excitation of the harmon- 



ics, see below. We, therefore, develop a kinetic theory of 
a magnetized plasma at finite temperature which incor- 
porates strong correlations. 

Theory. We start from the retarded longitudinal di- 
electric function e l of a one-component correlated plasma 
in a magnetic field [T9] . 

e\w,k) = ^e l3 (oj,k) = l-V(k)Il l (uj,k), (2) 

where e^- is the dielectric tensor, and V(k) = 2ttQ 2 /{k 2 + 
ft 2 ) 1 / 2 is the Fourier transform of the 2D Yukawa poten- 
tial (our results include the Coulomb case, for n — >• 0). 
The longitudinal polarization Ii l contains all interaction 
effects and can be expressed through the uncorrelated 
(Vlasov or RPA) polarization IIo via the Bethe-Salpeter 
equation [20] or via 

rf(a,,fc;B) = — — 3^§^T m , (3) 
1 + V (^)ILq(cl;, k; B)G{k, uj; B) 

where G(k,u;B) is the so-called local field correc- 
tion [21. However, G is known only approximately, re- 
sults for the electron gas have been given e.g. by Singwi 
et al. [21 j. For a strongly correlated plasma Kalman and 
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FIG. 3: (Color) MD collective excitation spectra, L(k,u) + T(/c,cj) (color) vs. theoretical dispersions, Eqs. Q, based on 
the corrected polarization fto (black and white lines). Left column: F — 200, Right column: F — 150 (e, f) and F — 600 (g, h). 
Note that frequencies are scaled in units of cj c . 



Golden observed [10] that the collective modes in the cold 
fluid description are well reproduced by a static approx- 
imation of Q in which the magnetic field is neglected, 
G{k,uj\B) « —D(k)/ujp(k), where D is the dynamical 
matrix of QLCA (a functional of the pair distribution 
function) [9 j. Here we extend this concept to finite tem- 
peratures and to the Bernstein modes (BM). To this end 
we compute the Vlasov polarization IIq (a;, k\ B) of a 2D 
Maxwellian plasma with Yukawa interaction in a mag- 
netic field, which is analogous to the quantum RPA re- 
sult US]: 

n&(«, k; B) = -?^e- £ „ I n (z), (4) 

where I n denotes the Bessel function, z = (kvT/u c ) 2 , and 
Vji = ksT '/m is the square of the thermal velocity. Using 
this result and Eqs. ([2| and Q the dispersion relation of 
the longitudinal modes is obtained from c l (u, fc; B) — 0, 

— z 00 2 

n=l v c/ 

where &p(k) = oo 2 (k) + D(k), and we use the approxi- 
mation D(k) = D L (k) + D T (k). Eq. ^ can be further 



simplified by noticing that z = (ka) 2 /(2f3 2 F) <C 1 for 
the situations of interest. For example, for /3 = 1 and 
T = 200, as in Fig. [3} z < 0.16, up to ka = 8. Then, us- 
ing I n (z) « (z/2) n /n!, Eq. ([5| yields the coupled MP and 
BM, which are modified by correlations. Approximately 
we find, for the MP, the well-known dispersion 

w?(fe) « c c 2 +^( fc ) ) (6) 

and, additionally, for the BM (n > 2), the frequencies 

o 1 / 7 \ n— 1 

col(k) « (n Wc ) 2 + [n^ p (fc)] 2 - (-) . (7) 

The latter are found to be undamped and to exist in the 
entire wave number range with frequencies rather close to 
the simulation results (note that since z <C 1 the second 
contribution to uj n is negligible for the shown parame- 
ters), cf. Fig. [3]i,h. This clearly confirms our interpreta- 
tion in terms of Bernstein modes. 

However, a closer inspection of the simulation data 
reveals three qualitative differences: 1) All BM are 
damped. 2) They exist only beyond a finite wave number 
kcr^ which increases with n, 3) The BM are not spaced 
by uj c but by a different frequency, cf. Fig. (3}},f. In fact 
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we find that the frequencies are very well described by 

(*l(k) « n 2 ^ uj 2 ^ = uj 2 + 2o;|, (8) 

where cji )00 is the large fc asymptotic of the MP ([6| and 
cje is the Einstein frequency describing the average os- 
cillation frequency of a test particle in the frozen envi- 
ronment of the other particles. The surprising absence of 
any wave number dispersion in the numerical results, cf. 
Figs. [2] and [3J suggests that BM excitation is governed 
by processes on the smallest length scale of single par- 
ticles. In fact, in a strongly correlated plasma particles 
are trapped in instantaneous local potential minima cre- 
ated by all other particles [22 j. The trajectories of these 
randomly oscillating particles are "bent" by the Lorentz 
force onto cyclotron orbits around the magnetic field di- 
rection giving rise to an average gyration frequency u)\ )00 . 
In the presence of thermal fluctuations nonlinear energy 
exchange with short wavelength magnetoplasmons leads 
to particle gyration with frequencies ncji 500 , connected 
with plasmon emission of these harmonics. Thus, the 
observed oscillations are "dressed Bernstein modes" aris- 
ing from a combined effect of magnetic field and strong 
correlations which increases the fundamental frequency 
from uo c to cji j00 . 

It is straightforward to present a corrected polariza- 
tion IL l ((jj,k]B) which reproduces all the above effects. 
The main difference, compared to Eq. Q, is that the de- 
nominators of the terms n > 2 are of the form [23] Cjuj — 
(ncji )00 ) 2 , which yields the correct BM dispersion ([8|. 
Further, by introducing a phenomenological damping 
constant v > 0, with Co = uu + iv, it is possible to correctly 
reproduce the damping of the BM and the finite k- values 
at which they emerge. This damping emulates collisional 
dissipation effects, which are missing in the Vlasov the- 
ory and in the QLCA, but are present in the simulations. 
Replacing, in the dispersion relation |H]), U l (u , k] B) by 
IL 1 {lj , k] B) yields, for the critical wave number of the 

n-th BM, 4 n ' fl = 2/^r 1 / 2 [(n- l)\uJ hoo i//w|] 1/(2n " 2) . 
Using v as a free parameter allows to reproduce the fi- 
nite values kffi and their increase with n observed in 
the simulations, as is shown by the dashed white lines in 
Fig. [3J We underline that this agreement is representa- 
tive for magnetized 2DYP in the strongly coupled liquid 
and crystal phases, what is supported by additional re- 
sults for different screenings, see right column of Fig. [3j 
In summary, we have demonstrated that a strongly cor- 
related 2DYP possesses, in addition to the well-known 
magnetoplasmon and magnetoshear modes, a sequence 
of Bernstein- type modes. In contrast to the familiar sit- 
uation of weakly nonideal high-temperature classical and 
low-temperature quantum plasmas where the BM appear 
at multiples of cj c , in a correlated plasma the modes are 
renormalized by correlations and appear at harmonics of 
the correlated MP frequency (uj 2 + 2a;!;) 1 / 2 . While har- 
monics of the plasma frequency may be excited in a cor- 



related 2DYP even without a magnetic field [24], there 
they are a very strongly damped. It is the combined 
effect of correlations and magnetic field which yields the 
strong harmonics signal with the frequencies reported 
in this Letter. The predicted dressed Bernstein modes 
should be observable, for example, in trapped ions, non- 
ideal electron-hole plasmas in quantum wells and in dusty 
plasmas. Using e.g. dust particles of 0.5/xra diameter and 
a magnetic field of 20T gives rise to a strongly coupled 
weakly magnetized liquid state with T ~ 100 and ~ 0.1 
which should allow for resonant detection of mode n = 2. 
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